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Spectral diffusion arising from ^®Si nuclear spin flip-flops, known to be a primary source of electron 
spin decoherence in silicon, is also predicted to limit the coherence times of neutral donor nuclear 
spins in silicon. Here, the impact of this mechanism on nuclear spin coherence is measured as 
a function of ^®Si concentration using X-band pulsed electron nuclear double resonance (ENDOR). 

The nuclear spin echo decays show that decoherence is controlled by ^®Si flip-flops resulting in 
both fast (exponential) and slow (non-exponential) spectral diffusion processes. The decay times 
span a range from 100 ms in crystals containing 50% ^®Si to 3 s in crystals containing 1% ^®Si. 

These nuclear spin echo decay times for neutral donors are orders of magnitude longer than those 
reported for ionized donors in natural silicon. The electron spin of the neutral donors ‘protects’ the 
donor nuclear spins by suppressing ^®Si flip-flops within a ‘frozen core’, as a result of the detuning 
of the ^®Si spins caused by their hyperfine coupling to the electron spin. 


Donors in silicon have been considered for use in qu an- 
tum information since the early days of the field.Ef^ 
Donors have both electron and nuclear spins which can 
be manipulated independently, and both have been con¬ 
sidered for use as potential quantum bits (qubits). While 
donor electron spins have received a majority of the 
attention,l2Hllthe nuclear spins are capable of much longer 
coherence times.^^^ This characteristic was utilized in 
the original Kane proposal for quantum computinj^ and 
gained attention later for building a quantum memory.l^ 
While exceptionally long T 2 times of donor nucl ear spins 
in silicon have already been demonstrated,^^ the me¬ 
chanics of nuclear spin decoherence are not yet fully un¬ 
derstood. In this study, we focus on neutral donor 
nuclear spin decoherence arising from interactions with 
^®Si nuclear spins in the silicon host environment. 

Spectral diffusion due to spin 1/2 ^®Si nuclei is a ma¬ 
jor source of decoherence for donor electron spins in 

silicoiPHISl 

and is predicted to be a major source of deco¬ 
herence for donor nuclear spins as well.^ While the pre¬ 
dicted coherence time for neutral donor nuclear spins 
in natural silicon (containing 4.7% ^®Si) was 0.5 s, several 
experimental works reported much shorter times (from 
hundreds of microseconds to tens of milliseconds). 
Coherence times presented here and by Wolfowicz et 
show that the limit from ^®Si spectral diffusion is actually 
longer than inferred from those previous experiments. To 
resolve the role of ^®Si spectral diffusion, we measure neu¬ 
tral ^^P nuclear spin coherence times in silicon crystals 
with ^®Si concentrations ranging from 1% to 50%.^ 

We And an inverse linear dependence of ^^P nuclear 
spin coherence time on ^®Si concentration (/), ranging 
from 100 ms at 50% ^®Si to 3 s at 1% ^®Si. The nuclear 
spin coherence time is about 1 second in natural silicon 
at 1.7 K; close to predictions of central spin stochastic 
models.^ However, contrary to the predictions, the ob¬ 
served spin echo decays are non-exponential. The decay 


times are two orders of magnitude longer than those mea¬ 
sured for ionized donors in natural silicorfi^I^orin NMR 
experiments on degenerately doped silicon.E3IiIl Appar¬ 
ently, the electron bound to a neutral donor protects the 
nuclear spin coherence from ^®Si flip-flops W detuning 
nearby ^®Si nuclear spins (a ‘frozen core’ ).l32H2i] This pro¬ 
tection might not be required in high purity isotopically 
enriched silicon, with a low content of ^®Si, where ^®Si- 
induced spectr al diff usion is no longer a dominant source 
of decoherencel ^ l ^l^^ 

Four phosphorus-doped silicon crystals with different 
concentrations of ^®Si isotopes were used in this work 
(Table |T]). In all crystals the donor concentration was 
about 10^®/cm^ which is low enough to ensure that other 
decoherence effects arising from dipolar interactions with 
donor electron spins are small compared to the measured 
^®Si spectral diffusion effects. The pulsed ENDOR ex¬ 
periments were conducted using a Bruker Elexsys E580 
spectrometer. Nuclear spin coherence times were mea¬ 
sured using an electron-mediated nuclear spin Hahn echo 
experiment.^ The combination of microwave and rf pulses 
enable a superposition state to be created on the donor 
electron, transferred to the nucleus, manipulated on 
the nucleus, and then transferred back to the electron for 
readout. For temperatures below 5 K, when the electron 
Ti relaxation was longer than 10 s, a light emitting diode 
(LED, 1050 nm) was flashed for 20 ms after each pulsed 
experiment in order to accelerate electron spin thermal- 
ization between repeated measurements. The “tidy” rf 
pulse to achieve nuclear spin thermaliz ation was not re¬ 
quired in these nuclear T 2 experiments.l^*23| Most of the 

data shown were measured with a static magnetic held 
(^ 0.35 T) oriented along a [001] crystal axis. Other held 
orientations were also examined to test the orientation 
dependence of the nuclear spin coherence times. 

The Hahn echo decay for phosphorus donor nuclear 
spins in natural silicon (/ = 4.7%) at 1.7 K is shown in 
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Sample 

"®Si (%) 

"®Si (%) 

®°Si (%) 

®^P/cm® 


98.1 

1.2 

0.7 

0.67 X 10^® 

''®Si-5% 

92.2 

4.7 

3.1 

0.8 X 10^® 

^®Si-10% 

87.2 

10.3 

2.5 

2.9 X 10^® 

2^81-50% 

50.2 

47.9 

1.9 

1.2 X 10^® 


TABLE I. Four ®^P-doped silicon samples used in this 
work. In each sample the ^®Si concentration was deter¬ 
mined by secondary ion mass spectroscopy (SIMS), and 
donor concentration was determined from ESR spin count¬ 
ing and independently confirmed by instantaneous diffusion 
slope measurements.^ All crystals were float-zone with the 
exception of ^®Si-5% (natural Si) which was Czochralski. All 
crystals had a volume on the order of a few cubic millimeters. 


Fiff. OTa). This decay is non-exponential and can be best 
fit 


v{t) = exp 




( 1 ) 


where t is the time interval between 7r/2 and tt pulses in 
a Hahn echo experiment. This functional form contains 
two decoherence terms. T 2 can be associated with vari¬ 
ous decoherence processes, including Ti-related processes 
and a broad variety of spectral diffusion mechanisms in a 
fast-motional regime, while Tsd is associated wit h spec- 
tral diffusion processes in a slow-motional regime 
As we will discuss, fast and slow motional regimes in our 
experiments are defined by how the rates of ^®Si nuclear 
spin flip-flops compare to the overall rate of decoherence. 
The stretch parameter n is in the range between 2 and 

o llOI24l25lTfl 

The temperature dependence of the extracted nuclear 
spin T 2 and Tsd for phosphorous donors in natural sili¬ 
con is shown in Fig. Qb). The Tsd term could only be 
extracted below 5 K because the linear T 2 term domi¬ 
nated the decays at higher temperatures. As seen from 
Fig-Bb), electron Ti controls the nuclear T 2 at temper¬ 
atures higher than 6 K, following the expected limit of 
T 2 n ~ 2Tie within experimental errors.^ However, elec¬ 
tron Ti continues growing below that temperature, while 
both T 2 and Tsd saturate at around 1 s showing lit¬ 
tle temperature dependence down to 1.7 K. This weak 
temperature dependence is consistent with Si-induced 
spectral diffusion being a dominant decoherence process 
for nuclear spin below 5 

The dependence of T 2 and Tsd on ^®Si concentration 
provides further evidence that ^®Si flip-flops are a major 
source of decoherence in our samples. Nuclear spin echo 
decaw for all four samples from Table |l] are shown in 
Fig. ^(electron spin echo decay for donors in natural 
silicon is also shown for comparison). The extracted T 2 
and Tsd at 1.7 K are plotted against ^®Si concentration 
in Fig. showing a relatively inverse linear depen¬ 
dence for both times. Within the experimental errors, 
the parameter n, shown in Fig. [^b), stays constant at 
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FIG. 1. (color online) (a) Nuclear spin Hahn echo decays 
for neutral donors in natural silicon (/ = 4.7%) at 1.7 K 
with magnetic field (~ 0.35 T) oriented along [001]. (b) Tem¬ 
perature dependences of nuclear spin T2 and Tsd (black 
squares and red circles, respectively) and electron spin Tie 
times (blue diamonds) for neutral phosphorus donors in nat¬ 
ural silicon. 


about 2.5 for all concentrations from 1% to 50%. T 2 and 
Tsd were also measured with the magnetic field oriented 
at different angles with respect to the crystal axis. How¬ 
ever, no noticeable orientation dependence was observed 
within experimental errors (10%) in either T 2 or Tsd- 
Three other decoherence mechanisms must be consid¬ 
ered here as potential contributors to nuclear spin deco¬ 
herence at low temperatures. These processes have been 
found to be significant in decohering electron spins of 
neutral donors. All three mechanisms are related to dipo¬ 
lar interactions with electron spins of other donors. The 
first two processes are cases of spectral diffusion aris¬ 
ing (1) from Ti-driven flips of electron spins of nearby 
donors,l^^^ and (2) from electron spin flip-flops in nearby 
donor pairs.l^^^Both cases are much less effective in deco¬ 
hering nuclear spins than electron spins since their effect 
scales proportionally with nuclear and electron gyromag- 
netic ratios (^1/1,600 in the case of nuclei). Using 
the electron T 2 times reported in Ref. |3], assuming donor 
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FIG. 2. (color online) nuclear spin Hahn echo decays 
for phosphorus donors in silicon with different ^®Si concentra¬ 
tions, measured at 1.7 K and magnetic field along [001]. The 
^®Si concentrations (/) are indicated for each curve. Electron 
spin Hahn echo decay for phosphorus donors in natural silicon 
is shown for comparison. 

densities of 10^®/cm^ (as used here) and considering tem¬ 
peratures below 4.8 K, we can then estimate the nuclear 
T 2 and Tsd from spectral diffusion processes (1) and (2) 
to be longer than 200 s. Thus, processes (1) and (2) are 
too slow to explain our T 2 data in Fig. Ufa). 

The third dipolar-related process to be considered is a 
“direct” flip-flop process.^ This involves a spin flip-flop 
between an electron of a “central” donor and an elec¬ 
tron of a neighboring donor (this is in contrast to “in¬ 
direct” flip-flops described in (2) above). Unlike other 
dipolar-related mechanisms, the effect of direct flip-flops 
does not scale with gyromagnectic ratio, decohering nu¬ 
clear spins as fast as electron spins. Direct flip-flops have 
been reported to limit electron spin coherence to 0.6 s for 
donors at 10^'^/cm^ in isotopically-purifled ^®Si crystals 
(45 ppm of ^®Si).® However, the inhomogeneous broad¬ 
ening in our samples (Table [l]) is 100-900 /iT which is 
20-200 times broader than the 5 /rT found in the afore¬ 
mentioned 45 ppm crystalsTaking into account the 
donor density in our samples we estimate that the direct 
flip-flop contribution to nuclear decoherence is about 
3 s in our ^®Si-l% sample and is even longer (>10 s) in 
the other three samples. 

The effect of ^®Si spectral diffusion on nuclear spin 
coherence of neutral donors in silicon has been ex¬ 
amined theoretically in the framework of a central spin 
problem while modeliM ^®Si spin flip-flops as a classi¬ 
cal stochastic process.!^ For natural silicon a ^^P nu¬ 
clear spin coherence time of 0.5 s was predicted with 
the held oriented along [001]. This prediction is very 
close to what was measured at that orientation in this 
work. Simulations of nuclear spin coherence showed an 
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FIG. 3. (color online) concentration dependence of (a) 
spectral diffusion times T 2 and Tsd, and (b) spectral diffusion 
exponential power, n, for nuclear spins of phosphorus 
donors in silicon at 1.7 K. Some error bars in (a) are smaller 
than their symbols. 


approximately inverse dependence on /, which also cor¬ 
relates with our results. However, the theory predicted 
exponential T 2 decays with n = 1, while our experiment 
shows non-exponential decays. The predicted orientation 
dependence was also not observed in our experiment. 

It is instructive to compare the decoherence of ^^P 
donor electron and nuclear spins caused by ^®Si flip-flops. 
The electron and nuclear spin echo decays differ in two 
ways: (i) the nuclear spin echo decays are over 3 or¬ 
ders of magnitude longer than electron spin echo decays, 
and (ii) the nuclear spin echo decays contain both ex¬ 
ponential and non-exponential components, unlike the 
electron spin echo decays that are dominated by only the 
non-exponential term.E^ Both electrons and ^^P nuclei 
see the same bath of ^®Si nuclear spins, with the bath’s 
dynamics suppressed in the “frozen core” where the nu¬ 
clear spins are detuned by the donor electron .EsHUl xhe 
main difference between electron and ^^P nuclear spins 
is the strength of their interactions to the ^®Si spin bath. 
Contact hyperfine interactions for an electron spin are 
much stronger than dipolar interactions for a ^^P nuclear 
spin, therefore the same ^®Si bath decoheres the electron 
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spin faster than the nuclear spin. This difference in the 
coherence timescale explains (i) and is also the key to 
understanding (ii). 

^®Si flip-flops can cause fast or slow motional effects 
depending on whether the flip-flop rate is fast or slow 
compared to th e strength of the pair’s interaction with 
the central gpin.ESIllIIZIMl Equivalently, the fast and slow 
regimes can be descriminated by comparing the flip-flop 
period to the overall coherence timescale (2r) of the cen¬ 
tral spin.l^ There is a broad distribution of ^®Si flip-flop 
rates, with the fastest rates being ^ 100 Hz and 10 Hz in 
nearest and next-nearest neighbor pairs, and much slower 
rates in more distant pairs. All these rates correspond to 
times much longer than 2t (~ 600 fis) when measuring 
electron spin echoes. In this case all flip-flops are in a slow 
motional regime, causing slow spectral diffusion with 
non-exponential echo de cays as seen in experimeniP^ and 
understood theoretically! ^° l ^^ ^ For nuclear spins, on the 
other hand, the timescale of the experiments lies within 
the broad distribution of ^®Si flip-flop times. Thus, there 
are fast and slow flip-flopping pairs that contribute to 
the decoherence, and consistently both exponential and 
non-exponential components are present in the decays. 

The nuclear spin coherence time in natural sili¬ 
con presented here is longer than measured earlier for 
ionized donors or donors in degenerately doped silicon. 
NMR measurements of nuclear spin decoherence in 
degenerately doped silicon h ave fo und times about two 
orders of magnitude shorter.lSnU Ionized donors mea¬ 


sured with EDMR had a coherence time of 18 ms,l^ and 
single donors measured with an SET had a coherence 
time of 60 ms.^ These measurements of ionized donors 
are in agreement with cluster correlation expansion sim¬ 
ulations by Witzel et al. (^30 ms).I^The longer coher¬ 
ence time for our isolated neutral ^^P donors supports 
the “frozen core” pictur^i^lt^ where most ^®Si pairs near 
a central spin are too detuned by the donor electron spin 
to flip-flop. 

In conclusion, we have experimentally studied the ef¬ 
fect of environmental ^®Si nuclear spins on neutral donor 
nuclear spins decoherence in silicon. Two contributors 
have been resolved arising from fast and slow flip-flopping 
^®Si nuclear spin pairs. We find that both contribu¬ 
tions exhibit a linear dependence on ^®Si concentration. 
Our results demonstrate long coherence times for neu¬ 
tral donor nuclear spins, ranging from 100 ms in crystals 
containing 50% ^®Si to 3 s in crystals containing 1% ^®Si, 
and are in agreement with the picture that an electron 
bound to a donor protects the donor nuclear spins from 
^®Si flip-flops. 
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